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ABSTRACT 

Emission-line regions in active galactic nuclei and other photoionized nebulae 
should become larger in size when the ionizing luminosity increases. This "breath- 
ing" effect is observed for the H/3 emission in NGC 5548 by using H/3 and optical 
continuum lightcurves from the 13-year 1989-2001 AGN Watch monitoring campaign. 
To model the breathing, we use two methods to fit the observed lightcurves in de- 
tail: (i) parameterized models and, (ii) the MEMECH0 reverberation mapping code. Our 
models assume that optical continuum variations track the ionizing radiation, and 
that the H/3 variations respond with time delays r due to light travel time. By fitting 
the data using a delay map ^(r, F c ) that is allowed to change with continuum flux F c , 
we find that the strength of the H/3 response decreases and the time delay increases 
with ionizing luminosity. The parameterized breathing models allow the time delay 
and the H/3 flux to depend on the continuum flux so that, r oc -Fjf and Fjj/3 oc F". 
Our fits give 0.1 < < 0.46 and 0.5 7 < a < 0.66. a is cons i stent with previous 
work bv lGilbert fc Peterson! l)2003j) and lGoad. Korista fc Kniggd ll2004D Although we 
find /3 to be flatter than previously determined bv IPetersonet^th 1)2002]) using cross- 
co rrelation methods , it is closer to the predicted values from recent theoretical work 
bv lKorista fc Goadl 12001 . 

Key words: galaxies: active - galaxies: individual (NGC 5548) - galaxies: nuclei - 
galaxies: Seyfert 



1 INTRODUCTION 

Photoionization models predict the sizes of Stromgren zones 
for H II regions and planetary nebulae ionized by hot stars 
of various luminosities and spectral types. Higher luminos- 
ity can maintain a larger mass of ionized gas. Dynamical 
tests of photoionization models are rare. The ionizing stars 
evolve in luminosity at rates too slow for humans to directly 
observe changes in radius. However, active galacic nuclei 
(AGN) vary on much shorter time-scales (days to months). 
Rapid variations in the ionizing luminosity emerging from 
an AGN should cause the photoionized region to expand and 
contract. This 'breathing' of the emission-line region is an 
interesting test of photoionization models. 

Although the time-scales of the variations are conve- 
nient to human observers, unfortunately, the angular sizes 
of the broad emission-line regions are too small to be re- 
solved directly. A nebula 100 light days across at a distance 
of 100 Mpc spans only 180 micro-arcseconds. Fortunately, 
light travel times within the nebula introduce time delays 
for any changes in the line emission. The ionizing radia- 
tion from the innermost regions of an AGN is reprocessed 
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by gas in the surrounding broad line region (BLR). As the 
central source varies, spherical waves of heating and ion- 
ization, cooling and recombination, expand at the speed of 
light through the BLR. A change in ionization causes a cor- 
responding change in the reprocessed emission. In AGNs, 
the recombination timescales of the gas in the BLR is very 
short compared to the light travel time, so the delay seen 
by a distant observer is dominated by the light travel time. 
Hence, we see line emission correlated with the continuum 
but with a time delay, r. A gas cloud 1 light day behind 
the ionizing source will be seen to brighten 2 days after the 
ionizing source flux rises. Thus, we can use light travel time 
delays to measure the size of the region that is responding 
to variations in the ionizing flux, where the reverberation 
radius is (R) ~ (r)c. Echo map ping, or reverberation map- 
ping jBlandford fc McKedll982lh aims to use this correlated 
variability to determine the kinematics and structure of the 
BLR, as well as the mass of the central supe rmassive black 
hole fe.g. |Petersonll993l : lPeterson et al 12004 and references 
therein). 

The nearby (z = 0.017) Seyfert 1 galaxy NGC 5548 has 
been intensively monitored in the optical range for 13 years 
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(198 9-2001) by the intern ational AGN Watch 1 consortium 
(e.g. iPeterson et"aHl2002l) . Those data spanning the source 
in a wide range of luminosity states are ideal for searching for 
this 'breathing' effect. In this paper, we investigate the lumi- 
nosity dependence of the H/3 emission and present two meth- 
ods of fitting the data, accounting for the 'breathing' using 
(i) parameteri zed models and (ii) the reverberation mappin g 
code MEMECHO faorne. Welsh fc PetersoJl99ltlHorneil994j) . 
Although previous work b v |P et*eTsori^^iL| 1^9^1120021) : 
lOilbert fc Peterson! (|200.^ : lOoad et all l)2004h has studied 
the luminosity dependence of the H/3 emission in NGC 5548, 
this study applies alternative techniques to characterise the 
'breathing'. In ^ we describe the echo mapping techinque 
and discuss the luminosity dependence of the emission-line 
lightcurve. In Sj3]we present the luminosity dependent pa- 
rameterized models followed by the MEMECHO method in 
The results of these methods of fitting the data, and their 
implications, are discussed in i|H]and we summarize our main 
findings in SjEl 



2 ECHO MAPPING 

The emission line flux, Fi(t), that we see at each time, 
t, is driven by the continuum variations, F c (t), and arises 
from a range of time delays, r. The emission line lightcurve 
is therefore a delayed and blurred version of the contin- 
uum lightcurve. In the usual linearized echo model the line 
lightcurve is modelled as 



F c (t) = F c + AF c (t) 



F t (t) = F,+ 



V (r) AF c (t ~ r)d,T 



(1) 
(2) 



where $(r) is the transfer function, or delay map. We can 
adopt a continuum background level F c , somewhat arbitrar- 
ily, at the median of the observed continuum fluxes. F] is 
then a constant background line flux that would be pro- 
duced if the continuum level were constant at F c . 

A simple way of determining the size of the emission- 
line region is to determine the time delay (or 'lag') be- 
tween the line and continuum lightcurves using cross- 
correlation. Taking the centroid of the cross-correlation 
function (CCF) as t he lag gives a luminosi ty-weighted 
radius for the BLR ijRobinson fc Perezl ^)90). However, 
the cross-correlation function is a convolution of the de- 
lay map, vE^r), with the auto-correlation function (ACF) 
of the driving continuum lightcurve. It is therefore pos- 
sible that changes in the measured cross-correlation lag 
arise from changes in the continuum aut o-correlation func- 
tion r ather than in the delay map (e.g. |Rolpm5gii_fc_Fjerej 
119901: iPerez. Robinson fc de La Fuentdll992t IWelshlll999T 
If the continuum variations become slower, a sharp peak 
at low time-delay in the delay distribution will be blurred 
by the broader auto-correlation function and the peak of 
the cross-correlation funct ion will be shifted to larger de- 
lays dNetzer fc M aoz 199(j). A typical delay map may have a 
rapid rise to a peak at small lag, and a long tail to large lags. 
The asymmetric peak in 9(t), will shift toward its longer 
wing when blurred by the auto-correlation function. Thus 

1 http: / /www.astronomy.ohio-state.edu/'~agnwatch/ 



the lag measured by cross-correlation analysis depends not 
only on the delay distribution, ^ (r), but also on the char- 
acteristics (ACF) of the continuum variations. 

Previous analysis of t he AGN Watch data for NGC 5548 
by iPeterson et alJ <|2002T) determined the H/3 emission-line 
lag relative to the optical continuum, on a year by year ba- 
sis, using cross-correlation. These authors find that the lag 
increases with increasing mean continuum flux. To improve 
upon the CCF analysis we use the echo mapping technique 
to fit the lightcurves in detail. However, the linearized echo 
model (Eq. |5J is appropriate only when the delay map is 
independent of time (static). In this paper we extend the 
model to search for changes in the delay map with luminos- 
ity. 



2.1 Luminosity-Dependent Delay Map 

The above linearized echo model (Eq. |2Jl assumes that the 
line emissivity and continuum flux can be related by a lin- 
ear function and thus is appropriate only for responses that 
are static. In principle, the delay map may change with 
time, for example, due to motion, or changes in quantity, 
of line-emitting gas within the system. The delay map may 
also change with lumin osity. In the "local op timally emitting 
clouds" (LOC) model fealdwin et alJll995li at each time de- 
lay there is a variety of gas clouds with differing properties, 
and those most efficient at reprocessing tend to dominate 
the line flux emerging from the region. A change in ionizing 
luminosity induces a change in the efficiency of reprocessing 
at each place in the region and so the time delay at which 
the line emission is dominant will change. When a cloud is 
partially ionized its response may initally be large so that 
increasing luminosity increases the depth of the ionized zone 
on the face of the cloud. Once the cloud becomes completely 
ionized, however, further increases in ionizing flux are less ef- 
fectively reprocessed. The line flux saturates, and may even 
decrease with increasing ionizing flux due to either ioniza- 
tion or decline in the recombination coefficie nts caused by an 
incre ase in gas temperature llO'Brien. Goad fc Gondhalekarl 
Il995l) . 

To account for these effects, we generalise the echo 
model by allowing the delay map to be luminosity- 
dependent, ^(t, F c ). The response we see at time t from 
a parcel of emission-line gas located at time delay r is set 
by the luminosity of the nucleus that we saw at the earlier 
time t — t. Thus, 



Ft(t)=Fi + 



*[T,F c (t-r)] AF c (t-r)dT 



(3) 



2.2 Luminosity dependence of H/3 flux 

A well-established correlation between continu um and 
emiss i on-line properties is the 'Bal dwin' effect fealdwinl 
Il977t lOsmer. Porter fc Greer] Il994l) where, in different 
AGN, broad emission- line equivalent width is observed to 
decrease with increasing continuum level. The relationship 
between the line luminosity, Li, and the continuum luminos- 
ity, L c , can be described by 



Li oc 



(4) 
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find that a w 0.83 for 



iKinnev. Rivolo fc Koratkarl 
CIV, and a « 0.88 for Lya. 

Within a single source, various studies have shown that 
emission li nes have a nonlinear response to continuum vari- 

iDietrich fe Kollatschnv 



19951). This 'intrinsic Baldwin effec 


' jKinnev et al.lll990l 


Krolik et alJll99ll iPoeee & Petersonl 


ll99llKorista & Goadl 


2004) where the H/3 emission- line response to variations in 



the continuum decreases with in creasing continuum level 
has been observed for NGC 55 48 (iGilbert fc PetersodElOol 
iGoad. Korista fc Knigge1l2004h . In terms of the continuum 
flux, F c , and the line flux, Fi, the intrinsic Baldwin effect is 
described by 



Fi oc Fc 



(5) 



This nonlinearity can be seen by simply examining the 13- 
year lightcurves (Fig. 0. In the lowest state (1992), the 
trough in the H/3 lightcurve is deeper than in the contin- 
uum lightcurve, whereas in the highest state (1998), the H/3 
peak is not as pronounced as the continuum. This argument 
neglects that light travel time delay smears out the emission- 
line response, but this should happen to both the peak and 
the trough. 

The relation between the optical continuum flux at 
5100A and the H/3 line flux is examined in Fig.|5](a), where 
a power-law (Eq. |SJ with a = 0.63 gives a good fit. Here 
we corrected the optical continuum flux at 5100 A, F opt , 
for the background host galaxy contribution, F ga i, where 



we take F, 



3.4 x 10" 



A as deter- 
mined bv iRomanishin et alJ dl995l) . The narrow-line com- 
ponent of the H/3 line has also been removed, and we use 
FHB(narrow) = 6.7 x 10~ 14 erg s _1 cm -2 determined by 



IGilbert fc Petersonl i2003f) . We included a time delay of 
17.5 days between the continuum flux and H/3 flux to re- 
move reverberation effects which was determined by cross- 
correlation of the full 13-year lightcurves. This time delay 
was subtracted from the times of each of the H/3 data points 
and the continuum flux at this new time det ermined via 
linear interpolation. IGilbert fc Petersonl i2003|) do a more 
detailed analysis allowing for the different time delays ob- 
served each year, and when adopting the same galaxy contin- 
uum backgroun d and H/3 narrow-l ine component, determine 
a = 0.65±0.02. lGoad et al.l(l2004h find that the slope of this 
relation is not constant, but decreases as the continuum flux 
increases, an effect which is predi cted by the photoionization 
models of lKorista fc Goadl ll2004h . However, the driving ion- 
izing continuum maybe closer to that observed in the UV 
at 1350 A, and so previous observations of NGC 5548 at 
this wavelength by the AGN Watch using IUE and HST 
JClavel et alJll99ltlKorista et alJll995h can be used to cor- 
rect the relationship s determined by the op tical continuum. 
Using the IUE data Petersoiie^^Jjj2002) finds a relation 
Fopt oc F^ v 6 , while jGilbert fc Peterson! (I2003T) find a rela- 
tion of F op t oc F[ 
we find a relation F op t oc F' x 
no time delay between the optical and UV continuum and 
linearly interpolating to get the optical continuum at the re- 
quired times. From this result, we get a relationship between 
the H/3 flux and the ionizing UV flux of F H (j oc F^y 1 . How- 
ever, the UV flux can be combined directly with the Hf5 flux 
to determine this relationship. Including cross-correlation 
time delays for the relavent years (19.7 days for the 1989 



Combining both IUE and HST data, 
02 (see Fig. 121(b)) assuming 
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Figure 1. Lightcurves for optical continuum flux (10~ 15 erg 
s — 1 cm _2 A , at 5100 A) and H/3 emission line flux (10 — 13 erg 
s" 1 cm" 2 ) from the 1989-2001 AGN Watch data on NGC 5548. 
LS and HS marks the lowest (1992) and highest (1998) states 
of the lightcurves respectively. Sol id lines mark the galaxy con- 
tribution to the continuum flux (Romanishin ct al. 1995) and 
the n arrow- line contribution to the H/3 flux iGilbert fc Petersonl 
2003). Dashed lines separate the observing seasons. 



data and 13.6 days for the 1993 data), we determine this 
relation directly to be F H p oc F% v 6±om (see Fig.H(c)). 



2.3 Luminosity dependence of time delay 

As the ionizing luminosity varies we expect the size of the 
photoionized region to expand and contract - a larger lu- 
minosity should ionize gas to a greater distance. We now 
consider a couple of simple theoretical predictions for this 
effect. If the BLR acts as a simple Stromgren sphere with 
uniform gas density, then one would predict that R oc L 1 ' 3 . 
If instead we assume that the response in an emission line 
will be greatest at some density, n, and ionization pa- 
rameter, F oc L/R 2 n, then it is easy to show this pre- 
dicts R oc L 1 / 2 (for a particular value of the product Tn) 
iPeterson et al.lf2002f) . More de tailed photoionization mod- 
eling (using the LOC model) bv lKorista fc Goadl i2004Tl pre- 
dicts a responsivity-weighted rad ius scaling as R oc L 23 for 
H/3. Photoionization model s by lO'Brien et alJ lll995l) and 
also iKorista fc Goadl <2004h both come to the conclusion 
that a relationship between emission-line lag and incident 
con tinuum level is due to a non-linear emissio n-line response. 

[ Peterson et all il999f) and more recentlv |Peterson et alJ 
l|200^1 used a year by year cross-correlation analysis of 
the AGN Watch data for NGC 5548 to show that the H/3 
emission-line lag (relative to the optical continuum) is cor- 
related with the mean optical continuum flux (at 5100A). 
As the mean optical continuum flux increases, the lag, and 
hence the size of the H/3 emitting region, is seen to in- 
crease. Using t he ful l 13-year lightcurves for NGC 5548, 
IPeterson et alJ (l2002h find r oc Fo P t 5 , though with much 
scatter. They argue, however, that the UV continuum (at 
1350A) is much closer to the driving ionizing continuum 
than the optical continuum used. Correcting for the re- 
lationship between the optical and UV continuum (using 
Fopt oc Fuy 3 ) leads to r oc F^ v ° as predicted assuming that 
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Figure 2. (a) Narrow-line subtracted H/3 flux vs. optical continuum flux at 5100A (F op t), with the host galaxy continuum flux removed 
(F ga i ) for the 13-year data set. The solid line indicates Fjjp oc F® p ^ 3 . A time delay between the continuum and H/3 fluxes of 17.5 days 
is inc luded, (b) UV continuum flux at 1350A vs. F op t - F ga i. UV data are taken from the IUE campaigns in 1989 (circles: IClavel et alj 
1991) and 1993 (squares) as well as the HST campaign in 1993 (triangles: iKorista et al1ll99fJ) . The solid line indicates the best fit of 
F op t oc F^'^? 02 . (c) Narrow-line subtracted H/3 flux vs. UV continuum flux at 1350A (symbols as in (b)). A time delay between the 
H/3 and UV fluxes has been added with the H/3 fluxes from 1989 19.7 days behind the UV and for the 1993 data the delay is 13.6 days. 
The solid line indicates the best fit of F H/3 oc F^y 8±0 01 . 



the emission will be greatest at some particular gas density 
and ionization parameter. 

To test the predictions, using more complex methods 
than the cross-correlation, we have fitted the data allow- 
ing for these breathing effects. Firstly, we present our pa- 
rameterized models and then the MEMECHO fits to the 13- 
year (1989-2001) AGN Watch optical continuum and H/3 
lightcurves for NGC 5548 before discussing these results, 
their findings and implications. 



3 PARAMETERIZED MODELS 

In this method we model the delay map and include param- 
eters to allow it to be luminosity dependent. We choose to 
model the delay map as a Gaussian in lnr, 



*W = T exp "5 



ln(r/r ) 



Alni 



(6) 



where to is the peak of the Gaussian and A In r is the width 
of the Gaussian, scales the strength of the delay map. 
We include a normalisation factor, l/A, chosen so that $o = 
J ^ (r) dr, where, 



.4 



exp 



2tt (A In r 




(7) 



We select this model as it ensures causality ("if (r) = for 
t < 0) while allowing parameters to and A lnr to control 
the centroid and width of the delay distribution. 

We introduce three 'breathing' parameters to allow the 
delay map to be luminosity dependent. To account for the 
two effects discussed in H2.2I and >I2.3I we introduce a to 



model the 'intrinsic Baldwin effect', so that Fh@ oc F" pt , and 
to allow the mean time delay of the delay map to depend 



on the continuum flux i.e., r oc F^ pt . A further parameter, 7, 
is introduced to allow the width of the delay map, A In r, to 
depend on the continuum flux, so that A lnr oc t . Fig- El 
illustrates how these parameters affect the delay map as the 
continuum flux increases. A linear line response corresponds 
to a — 1.0, whereas /3 = 0.0 means that there is no depen- 
dence of time-delay with luminosity, and 7 = 0.0 ensures 
A lnr is constant. A positive value for indicates that the 
time delay is increasing with increasing continuum flux, and 
similarly, a positive value for 7 indicates that the width of 
the delay map is increasing with increasing continuum flux. 

One complication is the host galaxy's contribution, 
Fg a i, to the observed optical continuum flux F op t(t). We 
treat this by defining the normalised lightcurve (with the 
background galaxy contribution subtracted) as 



X(t) = ^ 



Fgai 



(8) 



where F ga i is the background galaxy continuum flux and 
Fopt is the mean continuum flux, so that X = 1 when F op t = 
Fopt, and X — when F op t = F aa i- In all fits we adopt a 
background galaxy continuum flux, F ga i = 3.4 x 10~ 15 erg 



A" 



determined bv lRomanishin et al.l (Il99 



We include the three 'breathing' parameters (a, and 
7) to our model as follows, 



*o = *o [X(t 
to = to [X(t- 



Alnr 



A lnr [X(t-T)Y> 



(9) 
(10) 

(11) 



where $0, to and A lnr are just the values of ^0, to and 
A lnr at the mean optical continuum flux, F op t- The H/3 
lightcurve is then just given by 

F H0 (t)=F^+ / 9[T,X{t-T))AFopt{t-T)dT,Q2) 
Jo 

where, 
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Figure 3. Delay maps for different values of the 'breathing' parameters a, (3 and 7. Different lines indicate different continuum fluxes 
where F opt = 5.0 (solid), 10.0 (dashed), 15.0 (dot-dashed) and 20.0 (dotted) xl0~ 15 erg s" 1 cm" 2 A -1 , tq = 17.5 days and Alnr = 0.3 
throughout. 



F op t(t) — F op t + A_F op t(t) 



(13) 



As the delay map is Gaussian in lnr, not r, the delay 
map is asymmetric and has a long tail to high time delays 
(unless Alnr <C 1). Although we have parameterized the 
delay map in terms of To, the lag at which the delay map 
peaks, it is useful to characterise the delay map in terms of 
the median lag, r me( j, and the mean lag, (r), as these values 
are more directly comparable to the cross-correlation lags. 
We define these quantities as 

1 _ f m ' d ^{r)dr 

2 / °°*(r)dr ' 

and 



(14) 



J °° rfr (r) dr 
JT*(r)<*r 



(15) 



We now present tests of a series of models allowing the 
delay map to vary in different ways with continuum flux. The 
models are detailed in sections l3.1l to l331 with Fig.|7|showing 
the luminosity-dependent delay maps recovered from these 
models. The parameters found from the fits are detailed in 
Table [5] Initially, the delay map is made to be static. 



3.1 Static (S) 

The importance of allowing the delay map to be luminosity- 
dependent is highlighted when fitting the data with a static 
delay map. In this static model the delay map is independent 
of continuum level, thus we fix a = 1.0, /3 = 0.0 and 7 = 0.0. 
The free parameters are 9o, Tb~, Alnr and Fjjg, which are 
adjusted to give the best fit, determined by minimizing x 2 - 
The uncertainties are derived by A\ 2 = Xmin/N- The delay 
map is convolved with the continuum lightcurve to give the 
predicted line lightcurve (see Eg. 1121 with r max = 100 days. 
However, we need to know the continuum lightcurve at all 
times, and so we linearly interpolate between the continuum 
data points for this purpose. In all the parameterized model 
fits to the 13-year data set we take F op t = 9.73 x 10~ 15 
erg s" 1 cm -2 A , the mean of the optical continuum flux 
from this data. Not suprisingly, this static model fits the 13- 
year H/3 lightcurve poorly, with x 2 /1248 = 9.6 (where 1248 
is the number of H/3 data points). Table [5] gives the best 
fitting values of the parameters for this model and the delay 
map is shown in Fig. [7] 



We examine how the delay map changes with contin- 
uum flux by fitting this static model to the data on a 
year by year basis (see Table . When fitting the static 
model to the separate years, it was found that there was 
often not enough information to clearly determine Alnr, 
therefore we fix Alnr = 0.66, the value that is deter- 
mined from the global fit to the 13-year data set. Fig. 2] 
shows the change in the delay map between the lowest state 
(Kp~ t = 6.7 x 10" 15 erg s _1 cm -2 A" 1 in 1992) and the 
highest state (F op t = 13.5 x 10~ 15 erg s _1 cm~ 2 A -1 in 
1998). The delay peak increases by a factor of ~ 3 from 5.7d 
to 18. Od in time delay, while the height of the peak drops 
a factor of ~ 5 from 0.1 to 0.02. Thus, the H/3 response 
(the area under the curve) decreases by a factor ~ |. In 
this figure the cross-correlation function (CCF) and contin- 
uum auto-correlation functions (ACF) are shown for both 
of these years for comparison with the delay map. The CCF 
peak is close to the median of the delay map. In calcula ting 
the CCF and ACF we used the lWhite fc Peterson! lll994t) im- 
plementation of the interpolation cross-correlation method 
llGaskell fc Sparkdll986l iGaskell fc Petersonlll987l) . 

Fig. |S] shows how the time delay increases with in- 
creasing continuum flux. We have plotted r me£ j as this is 
less biased by long asymmetric tail of the delay map. Fit- 
ting a power-law of the form r med oc F^ pt to thi s gives 
(5 = 0.92 ±0.16. This is similar to the lPeterson et alJ (120021) 
findings where the lag was determined from cross-correlation 
and the resulting slope, j3 = 0.95 (see their Fig. 3). Thus, we 
confirm that the time delay increases with mean continuum 
flux. 



3.2 Bl 

In this model we allow two of the 'breathing' parameters, a 
and j3 to be free, while still fixing 7 = 0.0. The delay map is 
now luminosity dependent - the strength of the H/3 response 
and the time delay can vary with continuum flux. The free 
parameters in this model are a,/3,\f r o, to, Alnr and Fh0- 
As in the static model, we linearly interpolate the continuum 
lightcurve to get the continuum flux at the required times. 
This model fits the data significantly better than the static 
model with x 2 /1248 = 6.0 (Table0. The delay map (Fig.0 
clearly shows an increase in the time delay and a decrease 
in the H/3 response with increasing continuum flux. We find 
a = 0.66 ± 0.03 and /3 = 0.28 ± 0.05 for this model. It is 
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Table 1. Parameters for fits of the Static model to the yearly H/3 lightcurves of NGC 5548 between 1989-2001 as well as the global fit 
to all 13- years of data. From the global fit we determine Alnr = 0.66. In the static fits we fix AlnT to this value. 



Year 


F pt 


TO 


T 


(r) 


FfJ ft 




X 2 /dof 


1989 


9.92 ± 1.26 


12.7 ±0.8 


19.6 ± 1.2 


24.3 ± 1.5 


8.56 ± 0.03 


0.78 ±0.08 


198/ 132 = 


= 1.5 


1990 


7.25 ± 1.00 


11.4 ±0.7 


17.7 ± 1.1 


22.0 ± 1.4 


5.69 ± 0.04 


1.14 ±0.10 


278/94 = 


3.0 


1991 


9.40 ±0.93 


11.2 ± 1.1 


17.2 ± 1.7 


21.4 ± 2.1 


7.22 ± 0.05 


1.16 ±0.12 


88/65 = 
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1993 


9.04 ±0.90 


7.0 ± 1.2 


10.8 ± 1.8 


13.4 ± 2.3 


7.78 ± 0.03 


0.41 ±0.09 


386/142 = 


= 2.7 


1994 


9.76 ± 1.10 


10.9 ± 1.0 


16.8 ± 1.6 


20.9 ±2.0 


7.68 ± 0.03 


0.88 ±0.09 


324/128 = 


= 2.5 


1995 


12.09 ± 1.00 


10.8 ± 1.9 


16.7 ±2.9 


20.7 ±3.6 


9.35 ± 0.04 


0.59 ±0.11 


234/78 = 


3.0 


1996 


10.56 ± 1.64 


9.4 ± 0.5 


14.5 ±0.8 


18.1 ± 1.0 


8.08 ± 0.04 


0.61 ±0.09 


592/144 = 


= 4.1 


1997 


8.12 ±0.91 


12.5 ± 0.6 


19.3 ±0.9 


23.9 ± 1.2 


7.06 ± 0.04 


1.01 ±0.10 


384/95 = 


4.0 


1998 


13.47 ± 1.45 


18.0 ± 1.2 


27.7 ± 1.9 


34.2 ±2.3 


10.09 ±0.03 


0.65 ±0.09 


242/119 = 


= 2.0 


1999 


11.83 ± 1.82 


16.6 ± 1.3 


25.6 ± 2.0 


31.7 ± 2.5 


9.21 ± 0.03 


0.55 ±0.08 


99/86 = 


1.2 


2000 


6.98 ± 1.20 


6.0 ± 1.9 


9.3 ± 2.9 


11.6 ±3.6 


5.94 ±0.07 


0.81 ±0.19 


292/37 = 


7.9 


2001 


7.03 ±0.86 


10.6 ± 2.0 


16.4 ±3.0 


20.4 ±3.8 


5.29 ± 0.09 


1.31 ±0.26 


668/45 = 


14.8 


Global 


9.73 ±2.44 


11.4 ±0.6 


17.6 ± 1.0 


21.9 ± 1.2 


7.67 ±0.02 


0.79 ±0.01 


12016/1248 


= 9.6 
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Figure 4. Static delay map fit to the low state (1992, upper panel) and the high state (1998, lower panel). The central plots show the 
delay maps, where the median and mean time delays are marked by dashed and dot-dashed lines respectively. The mean time delay of 
the delay map in the high state is clearly larger, and the total response (the area under the curve) is lower, than in the low state. In 
these fits Alnr = 0.66. On the lightcurves, the dashed lines indicate the mean continuum and H/3 fluxes for the 13- year data set. Also 
shown (right) is the cross-correlation function (solid line) and continuum auto-correlation function (dashed line) for these years. 



interesting to note that a is lower than the value determined 
by fitting the static delay map to the yearly data (Fig.Q>Jk 



3.3 B2 

In an attempt to improve on model Bl we allow the width 
of the Gaussian delay map to vary as a function of contin- 



uum flux, letting 7 be a free parameter in the fit. Again, the 
continuum lightcurve is linearly interpolated to get the con- 
tinuum flux at the required times. This extra parameter re- 
sulted in a slight improvement of the fit, with x 2 /1248 = 5.9. 
Again, the delay map (Fig. |7J clearly shows an increase in 
the time delay and a decrease in the H/3 response with in- 
creasing continuum flux, a = 0.57±0.02 and /3 = 0.41 ±0.08 
for this model. 7 = — 0.20±0.06 allowing there to be a wider 
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Continuum Flux 

Figure 5. T me£ j vs. mean optical continuum flux for static delay 
map fit to yearly data. r me£ j (see Eq. 1141 is the median of the 
delay map. Solid line shows best fitting power-law to the data with 
Tmrd fx -F^ 2 ^ ' 16 . Dashed l ine indicates host galaxy continuum 
flux jp,omanishin et al. 199^). 



range of delays at lower continuum flux than at higher con- 
tinuum flux. 



3.4 B3 

The residuals of the B2 fit (Fig. |SJ, exhibit slow trends 
(timescale ~1 year) that are not fit by our model. There 
is no obvious correlation of these slow line variations with 
continuum flux, suggesting that they are due to a process 
that is independent of the reverberation effects. These trends 
may indicate a violation of the assumption that the distribu- 
tion of the line-emitting gas in the BLR is constant over the 
timescale of the data. However, as the dynamical timescale 
for the H/3-emitting gas, with R ~ re ~ 20 light days and 
An ~ 5000 km s -1 , is ~ 3 years, changes in the gas distribu- 
tion may well occur over the timespan of the observations. 
To allow for this in model B3 we fit a spline (with 26 nodes) 
to the residuals. Physically, this allows the background line 
flux to evolve with time e.g. accounting for differ ent amounts 
of lin e-emitting gas in the system. However, iGoad et alJ 
(2004) find that the index of the intrinsic Baldwin effect, 
a changes on these sorts of timescales and therefore these 
residuals might instead be interpreted as a consequence of 
that effect. Allowing the line background flux to vary im- 
proves the fit significantly, yielding x 2 /1248 = 3.1. The pa- 
rameters of the fit are similar to the B2 model (see Table|Sj, 
with a = 0.58 ± 0.03, f3 = 0.46 ± 0.07 and 7 = -0.24 ± 0.06. 

Fig.|H|shows the continuum and H/3 lightcurves for 1991- 
1992 and 1998-1999 and demonstrates how well the Static, 
B2, and B3 models fit the H/3 lightcurve. Particularly clear 
in this figure is the failure of the static model (dot-dashed 
line) to fit the deepest trough and largest peak. 



3.5 B4 

In the models so far we have used linear interpolation to de- 
termine the continuum flux at all the required times. This, 



however, can lead to unphysical lightcurves in the gaps be- 
tween the data points. It also takes no account of the noise 
in the data, so a wider delay map could be being recovered 
in an attempt to blur out the jagged continuum lightcurve. 
The reverberation mapping code MEMECH0 that we also use to 
determine the luminosity dependent delay map (see §0J uses 
maximum entropy methods to fit both the continuum and 
H/3 lightcurves. In this model we make use of the MEMECH0 
continuum fit to determine the continuum flux between the 
data points. The line background flux is still allowed to 
vary in the fit. The extra degrees of freedom allowed in the 
MEMECH0 fit to the continuum has again improved the fit, 
with xV 124 8 = 1-2- The width of the delay map from this 
model is seen to be thinner than the previous models (Fig.0 
because the MEMECH0 continuum lightcurve is smoother than 
the linearly interpolated continuum lightcurve. This has af- 
fected the value of 7, which is is now positive (see Table B^. 
and allows the width of the delay map to increase with in- 
creasing flux, the opposite to what was seen from the pre- 
vious model fits. The physical interpretation of a negative 
value of 7 is not obvious, and negative values of 7 found 
in B2 and B3 could be an artefact of using a linearly in- 
terpolated continuum model - a wider delay map at low 
continuum fluxes may be needed to smear out the jagged 
linearly interpolated continuum model. However, in the B4 
model the time delay still increases with increasing contin- 
uum flux, though the relationship is flatter (/3 = 0.10±0.01). 
a remains simliar to previous fits. 



4 MEMECH0 FIT TO THE NGC 5548 AGN WATCH 
1989-2001 LIGHTCURVES 

The echo mappin g computer code MEMECH0 (see lHorne et alJ 
1991; Home 1994, for technical details) can allow luminosity- 
dependent delay maps to be made by means of a maximum- 
entropy deconvolution of Eq. |3] The maximum-entropy 
method finds the simplest positive image that fits the data, 
balancing simplicity, meas ured by entro py, and realism, mea- 
sured, in this case, by x 2 <Horndll994Tl . 

The optical continuum lightcurve, F op t(t), must thread 
through the measured continuum fluxes, and the predicted 
H/3 emission- line lightcurve, Fnpit), must similarly fit the 
measured line fluxes. In our fit to the data points we ad- 
just the line background flux Fh/3, the continuum variations 
F op t(t), and the delay map $>(t, F opt ), we take F op t as the 
median of the continuum data. The fit is required to have 
a reduced /N = 1 ± y/2/N, where N is the number of 
data points. We require this to hold for the line flux mea- 
surements, and also for the continuum flux measurements. 
The continuum is split into several flux levels (indicated in 
the lower panel of Fig.|5J, and delay maps corresponding to 
each level are determined. When computing the convolution 
(Eq.[3J we linearly interpolate between these continuum lev- 
els to find the delay map that applies to the continuum flux 
at time t — r. Such an approach allows for fully non-linear 
line re sponses. Further details of this can be found in lHornel 

GUI). 

Fig-OUshows the measured lightcurves and the MEMECH0 
fit. The greyscale (middle panel) shows the delay map at 
each time. Fig. 1101 shows the luminosity-dependent delay 
map ^(TjFopt) (greyscale) reconstructed from the observed 
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Figure 6. Data and residuals for B2 model fit to the 1989-2001 H/3 lightcurve. Residuals plotted are (data — model), and \ = (data 
model) I a. 



Table 2. Parameters for fits to 1989-2001 H/3 lightcurve of NGC 5548. In all models the mean optical continuum level F opt = 9.73 X 10 -15 



erg s 1 cm 



-2 A" 1 . 





X 2 /1248 


a 


/3 


7 


TO 


Alnr 


Fh/3 


*o 


s 


9.6 


1.0 


0.0 


0.0 


11.4 ±0.6 


0.66 ± 0.05 


7.67 ±0.02 


0.79 ±0.01 


Bl 


6.0 


0.66 ± 0.03 


0.28 ±0.05 


0.0 


12.9 ± 1.0 


0.70 ±0.07 


7.98 ±0.02 


0.77 ±0.05 


B2 


5.9 


0.57 ±0.02 


0.41 ±0.08 


-0.20 ±0.06 


13.9 ±0.8 


0.56 ±0.05 


7.98 ±0.02 


0.68 ±0.01 


B3 


3.1 


0.58 ± 0.03 


0.46 ± 0.07 


-0.24 ±0.06 


12.3 ±0.7 


0.58 ±0.04 


7.97 ±0.01 


0.68 ±0.01 


Bl 


1.2 


0.64 ±0.01 


0.10 ±0.01 


0.46 ± 0.05 


16.4 ±0.3 


0.32 ±0.02 


8.05 ±0.01 


0.68 ±0.01 



lightcurves at each continuum level. The crosses indicate 
the median of the delay map and the dots indicate the up- 
per and lower quartiles (see Tab. The lower left panel 
projects the delay map along the time delay axis and thus 
indicates that the amplitude of the line response decreases 
with rising continuum level. The upper right panel gives the 
luminosity-averaged one-dimensional transfer function. At 
each continuum luminosity the range of delays is relatively 
narrow compared to the one- dimensional transfer function. 



As the luminosity rises, the mean delay increases. At mini- 
mum light the median delay is ~ 13 days, and this rises to 
~ 23 days at maximum light. Fitting a power-law of the form 
t ex F^ pt (with the galaxy background continuum removed) 
to the median time delay leads to f3 = 0.24 ± 0.08. 

In order to understand how robust the MEMECH0 re- 
covered luminosity-dependent delay map is we ran a short 
Monte-Carlo simulation. We generated ten sets of contin- 
uum and H/3 lightcurves with the data points shifted by a 
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Figure 7. Luminosity-dependent delay maps for the 5 parameterized models as well as the MEMECH0 recovered delay map. The reduced 
X 2 value of each fit is indicated. For the parameterized models, the dashed line indicates the dependence of the peak of the delay map on 
the continuum flux. For the MEMECHO model the dashed line indicates the best fit power-law, r oc F^ pt to the delay map. The dot-dashed 
lines indicate the minimum and mean continuum flux of the 13-year data, and the upper-limit of the plots is the maximum continuum 
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Figure 8. Continuum and H/3 lightcurves for 1991 - 1992 (left) and 1998 - 1999 (right). The linearly interpolated continuum flux model 
is marked with a solid line. The fit to the H/3 lightcurve can be seen for 3 of the different models (B3 (solid line), B2 (dotted line) and S 
(dot-dashed line)). 



Gaussian random number with a mean of zero and stan- 
dard deviation equal to the uncertainty on that data point. 
Ten different delay maps were then recovered using these 
lightcurves. The delay maps recovered are all close to the 
original. Fig. llll shows the ten realisations of the delay maps 
when the continuum flux is at the average of the low and 
high states (7.0 and 13.0 x 10~ 15 erg s~ x cm" 2 A -1 re- 



spectively). This shows the possible range of uncertainties 
expected in the delay map. 



5 DISCUSSION 

Both the parameterized fit and the MEMECHO fit to the 13- 
year continuum and H/3 lightcurves for NGC 5548 show that 
H/3 reverberations depend upon the continuum state in such 
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Table 3. Time delay at the median, lower and upper quartile of the MEMECHO luminosity-dependent delay map for each continuum flux. 



Continuum Flux 
(lCT 15 erg s" 1 cm" 2 A -1 ) 
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Figure 9. Lightcurves for optical continuum (bottom) and H/3 emission line flux (top) from the 1989-2001 AGN Watch data on 
NGC 5548. The MEMECHO fit is the smoothest model that achieves x 2 /N = 1 ± \J2/N for both lightcurves. In the middle panel, the 
greyscale shows the delay map at each time, the solid line marks the median of the delay map, where as the dotted lines mark the lower 
and upper quartiles. The horizontal dotted lines in the bottom panel show the flux levels at which each delay map is determined. 



a way that greater time delays occur for higher continuum 
states. From the parameterized breathing models we find 
that r oc t with j3 in the range 0.1 to 0.46, depending on 
the specific model, with the continuum at 5100 A corrected 
for the galaxy contribution. However, as noted in >I2.2I the 



driving ionizing continuum may be closer to that observed 
at 1350 A. Using our relationship that F opt oc F^ v 3 , the 
luminosity-dependent lag scales like r oc f^ 5-0,24 . The 
MEMECHO fit gives r oc F° p ? 4 , or r oc Fy v 3 with respect to 
the UV continuum. This would appear to be in contrast 
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Figure 11. 10 realisations of the MEMECHO recovered delay map 
at low and high continuum levels (7.0 and 13.0 X 10 — 15 erg s _1 
cm -2 A respectively) from Monte-Carlo simulations. 



with both the results of the Static model fits to the sep- 
arate years, where w e find r oc F® p ® 3 , and the findings of 
iPeterson et alJ (I2002T) where they find that CCF centroid 
lag scales as F^ 5 . Our values from the breathing models 
ar e closer to that p r edicte d by the photoionization models 
of iKorista fc Goadl i2004) who find that the responsivity 
weighted radius, R oc F®' 23 , although they noted that the 
power-law slope should differ somewhat in high (flatter) and 
low (steeper) continuum s tates. 

IPeterson et alJ (1200211 determine the lag on a yearly ba- 
sis using cross-correlation. As mentioned in § |2] the CCF is 
the convolution of the delay map with the ACF. Thus, it is 
clear that the time delay between the line and continuum 
light curves determined via cross-correlation depends on the 
shape of both the delay map and the ACF of the continuum. 
If the ACF is a delta function, then the CCF is identical to 
the delay map. Generally the ACF is broad and the delay 
map is asymmetric, thus the peak of the CCF will not nec- 
essar ily occur at the same time delay as in the delay map 
(e.g. IWelshlll999l) . The continuum variability properties of 
NGC 5548 vary from year to year leading to a change in 
ACF, and as an artifact of this, a change in lag could be 
measured without there being a change in the delay map. It 
is also the case that the accuracy of the CCF centroid de- 



pends on the length of the sampling window, the sampling 
rate, th e continuum variab ility characteristics and the data 
quality JPerez et al.lll992t) . Throughout the 13-year moni- 
toring campaign of NGC 5548 the mean sampling rate has 
varied, which could lead to a change in the determined lag. 
However, our results from the Static fits to the yearly data 
produce similar lags (when compar ing our T me d w i th th e 
centroid lag) and a similar slope to IPeterson et aD {2002), 
suggesting that the cross-correlation results may be measur- 
ing the true lags. 

Both the parameterized breathing fits and the MEMECHO 
fit to the full 13-year data show a much flatter dependence 
of the lag on the luminosity than the year-by-year analysis 
shows. There is scatter within the results from the yearly 
analysis. For example, there is a wide range of delay (~ 
9—18 days) when the continuum is in a low state with 
F op t ~ 7xl0 -15 erg s _1 cm~ 2 A , and there is also a wide 
range of continuum fluxes (7 - 12x 10 -15 erg s _1 cm -2 A ) 
that have a lag of approximately 17 days. A wide range of 
delays at low continuum flux and also the same lag at a 
wide range of fluxes favours a flatter relationship. It could 
be that these features dominate over the two years with the 
largest lag when fitting to the full 13-year data set. The 
tail to long delays on the delay maps determined by the 
parameterized models and the MEMECHO fits may also allow 
a flatter relationship. The difference in slope is certainly 
interesting and merits further investigation. 

From both the parameterized models and the MEMECHO 
fits, we also find that the amplitude of the H/3 response de- 
clines with increasing continuum luminosity, in other words, 
the change in H/3 flux relative to changes in the continuum 
is greater in lower continuum states . Deta iled photoioniza- 
tion calculations bvlKorista fc Goa dl j2004t ) predict this, and 
iGilbert fc Peterson! feOOSt and iGoad et alJ ll2004h find the 
same general result. From our parameterized fits to the data, 
we find Fh/3 oc F" pt wi th 0.57 < a < 0.66, which is consis- 
tent with that found bv lGilbert fc PetersonH2005t . Correct- 
ing our relationship to be relative to the UV continuum gives 
Flip oc F|° r ^ ~ 35 , which is c lose to the relationship shown in 
Fig 121 (c). IGoad et alJ <l2004l) find a range in this parameter 
(relative to the optical continuum) of 0.41 < a < 1.0, that 
correlates with the optic al continuum flux with a high degree 
of statistical confidence. IKorista fc Goadl (120041) find this to 
be a natural result, given that the H/3 equivalent width drops 
with increasing incident photon flux. Their model predicts 
slightly higher values for the responsivity to the ionizing flux 
(Fhp oc F® ^ 4 ~ ' 77 ) than we observe from the parameterized 
fits. 

The MEMECHO recovered delay map (Fig. Q panel (f ) and 
Fig. 1101 and the delay map from the B4 model (Fig. [T] panel 
(e)) seem similar, which is not suprising as they are both 
driven by the same continuum lightcurve model. However, 
the B4 model differs from the Bl, B2, and B3 models, partic- 
ularly in its value for /3, because of the differing continuum 
models used. MEMECHO allows the continuum model to have 
more freedom between the data points (compared to the 
linear interpolation used in models Bl - B3) resulting in a 
delay map with a smaller width at shorter time delays as 
less blurring of the continuum lightcurve is required. Nev- 
ertheless, the q parameter does not vary too much between 
the different models. It is also interesting to note that the 
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residuals in B4 are very flat, even without the spline fit, 
the long term trends seen in the residuals to B2 have dis- 
appeared. On close inspection of the MEMECHO fit to the con- 
tinuum (see Fig. 1 1 2> , it can be seen that between the data 
points the continuum model sometimes shows unphysical 
dips and peaks. For instance, if the model constantly dips 
down between the data points, the lag between the contin- 
uum and the line lightcurves can be shifted without altering 
the delay map. In this way it seems that MEMECHO accounts 
for the long timescale variations that cannot be described 
by the delay map in the B2 model. Although these spurious 
dips mean that the MEMECHO result may not be completely 
reliable, it is clear that the recovered lag luminosity rela- 
tionship is much flatter than found by the cross-correlation 
method and thus is an important independent check of the 
parameterized models. 

For the B3 and B4 models we allowed a spline fit to the 
residuals to account for slowly varying trends (timescale ~ 1 
year, see Fig.||jJ. We suggest that these long timescale vari- 
ations are not due to reverberation effects, but due to slow 
changes of the BLR gas distribution. In a study of the line 
profile variability of the first five years of the A GN Watch 
data on NGC 5548, IWanders fc Peterson! <1996ft find that 
although the H/3 emission line flux tracks the continuum 
flux, the H/3 emission line profile variations do not, and 
therefore are not reverberation effects within the BLR. In 
Fig. 1131 w e compare the lin e profile 'shapes' determined by 
Wan ders fc Peterson! il996l) for the fi rst five years of data 
with o ur residuals from the B2 model. IWanders fc Peterson! 
( 1996) define these emission-line 'shapes' to describe the 
relative prominence of features in the emission-line profile. 
Each profile is split into its red wing, core and blue wing and 
the 'shape' for each determined. A positive number for the 
shape indicates that it is more prominent than on average, 
whereas a negative number indicates that it is less promi- 
nent than on average. Although the shapes vary on similar 
timescales to the residuals, there is no obvious correlation 
with the profile shape variations. A comparison of the pro- 
file variability for the full 13-year data with our residuals 
warrants further study. 



6 CONCLUSION 

Our analysis of 13 years of optical spectrophotometric mon- 
itoring of the Seyfert 1 galaxy NGC 5548 from 1989 through 
2001 reveals that the size of the H/3 emission-line region in- 
creases as the continuum increases. Will also find that the 
strength of the H/3 response decreases with increasing con- 
tinuum flux. We have fit the H/3 emission-line lightcurve us- 
ing both parameterized models and the echo-mapping code 
MEMECHO, both methods show these effects. 

In our parameterized models, we allow the delay map to 
be luminosity-dependent. Our model is parameterized such 
that the H/3 emission can respond non-linearly to the op- 
tical continuum variations, i.e., Fnp oc F" pt . From our fits 
to the data we determine 0.57 < a < 0.66 which i s con - 
sist ent with previous fi ndings of lGilbert fc Peterson! J2003T) 
and lGoad et alJ <l2004h . However, the ionizing continuum is 
likely to be closer to the UV continuum (Fuv) than the 
optical continuum. Correcting for this (F opt oc Fy-y 3 ) gives 
Fuji oc F°y°~ - 35 . In addition, we allow the peak of the 



delay map, To, to be luminosity-dependent, to oc F% pt , and 
find 0.10 < f3 < 0.46. Correcting to be with respect to the 
UV continuum gives r oc i^° 5 ~ ' 24 . MEMECHO fits to the 
lightcurves also show these effects, and fitting a power-law 
of the form r oc t to the luminosity-dependent delay map 
gives t oc Fop? 4 , or r oc F^y 3 . The values we determine for 
(5 (corrected to be relative to the UV continuum) are not 
consistent with the simple Stromgren sphere (/3 = 1/3) and 
ionization parameter (/3 = 1/2) arguments and have a flat- 
ter slope t han the = 0.95 result from the cross-correlation 
analysis of iPeterson et al.1 (I2002I). Howeve r, they are close 
to the prediction of iKorista fc Goa d (2004) whose more de- 
tailed photoionization models predict /3 = 0.23. 

In our parameterized models we find slowly vary residu- 
als (timescale ~ 1 year) which we suggest are not reverber- 
ation effects, but indicate changes in the gas distribution. 
Comparison of these residuals with ve locity profile changes 
from the first five years of the data (IWanders fc Peterson! 
1996) are inconclusive. In future work we will examine ve- 
locity profile changes of the full 13-year lightcurves. 
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Figure 12. A section of the MEMECHQ fit to the continuum lightcurve covering the observing campaign in 1993 and 1994. 
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Figure 13. A comparison of the 'shapes' of the red wing, blue wing and core of the H/3 line profile as calculated bv lWanders &: Peterson! 
J1993) with the residuals of the B2 parameterized model, and the continuum and H/3 flux with the B2 model shown. 
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